Decoherence dynamics of a qubit coupled to a quantum two-level system 
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We study the decoherence dynamics of a qubit coupled to a quantum two-level system (TLS) in 
addition to its weak coupling to a background environment. We analyze the different regimes of 
behaviour that arise as the values of the different parameters are varied. We classify those regimes 
as two weak-coupling regimes, which differ by the relation between the qubit and TLS decoherence 
times, and a strong-coupling one. We also find analytic expressions describing the decoherence 
rates in the weak-coupling regimes, and we verify numerically that those expressions have a rather 
wide range of validity. Along with obtaining the above-mentioned results, we address the questions 
of qubit-TLS entanglement and the additivity of multiple TLS contributions. We also discuss the 
transition from weak to strong coupling as the parameters are varied, and we numerically determine 
the location of the boundary between the two regimes. 



I. INTRODUCTION 

There have been remarkable advances in the quest 
to build a superconductor-based quantum information 
processor in recent years 0,0,|ESlEIE0B0) EH 
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Coherent oscillations have been ob- 
served in systems of single qubits and two interacting 
qubits 0, H H 0. In order to achieve the desirable 
power of a functioning quantum computer, one would 
need to perform a large number of quantum gate oper- 
ations of at least hundreds of qubits. One of the main 
obstacles to achieving that goal, however, are the rela- 
tively short decoherence times in these macroscopic sys- 
tems (note that even a single superconductor-based qubit 
can be considered a macroscopic system). Therefore, 
there has been increasing experimental and theoretical 
activity aimed at understanding the sources and mech- 
anisms of decoherence of such s yste ms in recent years 

[3l3l3l3l3l3l3l3l3[3l2ai2a. 

The environment causing decoherence of the qubit is 
composed of a large number of microscopic elements. 
There is a large wealth of theoretical work on the so- 
called spin-boson model |22|. which models the environ- 
ment as a large set of harmonic oscillators, to describe 
the environment of a solid-state system. However, re- 
cent experimental results suggest the existence of quan- 
tum two-level syste ms ( TLSs) that are strongly coupled 
to the qubit [3 [3 III. Furthermore, it is well known 
that the qubit decoherence dynamics can depend on the 
exact nature of the noise causing that decoherence. For 
example, an environment composed of a large number of 
TLSs that are all weakly coupled to the qubit will gen- 
erally cause non-Markovian decoherence dynamics in the 
qubit (see, e.g., 0]). Note that the mechanism of qubit- 
TLS coupling depends on the physical nature of the qubit 
and TLS. The exact mechanisms are presently unknown. 

The effect on the qubit of an environment consisting 
of weakly coupled TLSs with short decoherence times is 
rather well understood. As was presented in Ref. [l8| . 
one takes the correlation functions of the TLS dynam- 



ics in the frequency domain, multiplies each one with a 
factor describing the qubit-TLS coupling strength, and 
adds up the contributions of all the TLSs to obtain the 
effective noise that is felt by the qubit. We shall refer 
to that approach as the traditional weak-coupling ap- 
proximation. In this paper we shall study the more gen- 
eral case where no a priori assumptions are made about 
the TLS parameters. We shall identify the criteria un- 
der which the traditional weak-coupling approximation is 
valid. We shall also derive more general expressions that 
have a wider range of validity, as will be discussed be- 
low. Furthermore, we study the criteria under which our 
weak-coupling results break down, and the TLS cannot 
be easily factored out of the problem. It is worth noting 
here that we shall not attempt to theoretically repro- 
duce the results of a given experiment. Although we find 
potentially measurable deviations from the predictions 
of previous work, we are mainly interested in answering 
some questions related to the currently incomplete un- 
derstanding of the effects of a TLS, or environment of 
TLSs, on the qubit decoherence dynamics. 

Since we shall consider in some detail the case of a 
weakly coupled TLS, and we shall use numerical calcula- 
tions as part of our analysis, one may ask why we do not 
simulate the decoherence dynamics of a qubit coupled to 
a large number of such TLSs. Alternatively, one may ask 
why we separate one particular TLS from the rest of the 
environment. Focussing on one TLS has the advantage 
that we can obtain analytic results describing the con- 
tribution of that TLS to the qubit decoherence. That 
analysis can be more helpful in building an intuitive un- 
derstanding of the effects of an environment composed of 
a large number of TLSs than a more sophisticated simu- 
lation of an environment composed of, say, twenty TLSs. 
The main purpose of using the numerical simulations in 
this work is to study the conditions of validity of our 
analytically obtained results. 

The present paper is organized as follows: in Sec. II 
we introduce the model system and the Hamiltonian that 
describes it. In Sec. Ill we describe the theoretical ap- 
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proach that we shall use in our analysis. In Sec. IV we use 
a perturbative calculation to derive analytic expressions 
for the relaxation and dephasing rates of the qubit in the 
weak-coupling regime and compare them with those of 
the traditional weak-coupling approximation. In Sec. V 
we numerically analyze the qubit decoherence dynamics 
in the different possible regimes. We also address a num- 
ber of questions related to the intuitive understanding of 
the problem, including those of qubit-TLS entanglement 
and the case of two TLSs. In Sec. VI we discuss the 
question of the boundary between the weak and strong 
coupling regimes, and we perform numerical calculations 
to determine the location of that boundary. We finally 
conclude our discussion in Sec. VII. 



II. MODEL SYSTEM 

We consider a qubit that is coupled to a quantum TLS. 
We take the qubit and the TLS to be coupled to their 
own (uncorrelated) environments that would cause de- 
coherence even if the qubit and TLS are not coupled to 
each other. We shall be interested in the corrections to 
the qubit decoherence dynamics induced by the TLS. The 
Hamiltonian of the system is given by: 

H = Hq + Htls + Hi + Hehv, (1) 

where H q and Htls are the qubit and TLS Hamiltoni- 
ans, respectively, Hi describes the coupling between the 
qubit and the TLS, and H Env describes all the degrees 
of freedom in the environment and their coupling to the 
qubit and the TLS. The qubit Hamiltonian is given by: 

£ q = -^i q) -^i q \ (2) 

where A q and e q are the adjustable control parameters 

of the qubit, and are the Pauli spin matrices of the 
qubit. For example, for the charge qubit in Ref. Q, A q 
and e q are the energy scales associated with tunnelling 
and charging, respectively. Similarly, the TLS Hamilto- 
nian is given by: 

tt A T LS ~ (TLS) £ TLS > (TLS) fn\ 

-Htls = — Y~ z ' ' 

The energy splitting between the two quantum states of 
each system, in the absence of coupling between them, is 
then given by: 

E a = ^A a 2 + eJ, (4) 

where the index a refers to cither qubit or TLS. For fu- 
ture purposes, let us also define the angles 8 a by the 
criterion 

tan# Q = (5) 



We take the interaction Hamiltonian between the qubit 
and the TLS to be of the form: 

Hl = -h^®a^\ (6) 

where A is the coupling strength between the qubit and 
the TLS. Note that the minus sign in Hi is simply a 
matter of convention, since A can be either positive or 
negative. It is worth mentioning here that the applica- 
bility of this form of interaction is not as limited as it 
might appear at first sight. Any interaction Hamiltonian 
that is a product of a qubit observable (i.e. any Hermi- 
tian 2x2 matrix) and a TLS observable can be recast in 
the above form with a simple basis transformation. 

We assume that all the terms in H^nv are small enough 
that its effect on the dynamics of the qubit+TLS sys- 
tem can be treated within the framework of the Marko- 
vian Bloch-Redfield master equation approach. It is well 
known that the effect of certain types of environments 
cannot be described using that approach, e.g. those con- 
taining 1// low-frequency noise. However, there remain 
a number of unanswered questions about the problem 
of a qubit experiencing 1// noise. For example, it was 
shown in Ref. 0] that the decoherence dynamics can 
depend on the specific physical model used to describe 
the environment. Therefore, treating that case would 
make it more difficult to extract results that are directly 
related to the phenomenon we are studying, namely the 
effect of a single quantum TLS on the qubit decoherence 
dynamics. We therefore do not consider that case. 

Depending on the physical nature of the system, the 
coupling of the qubit and the TLS to their environments 
is described by specific qubit and TLS operators. In prin- 
ciple one must use those particular operators in analyz- 
ing the problem at hand. However, since we shall present 
our results in terms of the background decoherence rates, 
which are defined as the relaxation and dephasing rates in 
the absence of qubit-TLS coupling, the choice of system- 
environment interaction operators should not affect any 
of our results. In fact, in our numerical calculations 
below we have used a number of different possibilities 
and verified that the results remain unchanged, provided 
that the background decoherence rates are kept constant. 
Furthermore, the background-noise power spectrum af- 
fects the results only through the background decoher- 
ence rates. Note that we shall not discuss explicitly the 
temperature dependence of the background decoherence 
rates. It should be kept in mind, however, that the back- 
ground dephasing rates generally have a strong temper- 
ature dependence in current experiments on supercon- 
ducting qubits. 

It is worth noting that, since we are considering a 
quantum-mechanical TLS, the model and the intermedi- 
ate algebra that we use are essentially identical to those 
used in some previous work studying two coupled qubits 
[^ESEii- However, as opposed to being a second qubit, 
a TLS is an uncontrollable and inaccessible part of the 
system. Therefore, in interpreting the results, we only 
consider quantities related to the qubit dynamics. 



III. THEORETICAL ANALYSIS: MASTER 
EQUATION 

As mentioned above, we take one particular element of 
the environment, namely the TLS, and do not make any 
a priori assumptions about its decoherence times or the 
strength of its coupling to the qubit. We assume that the 
coupling of the qubit to its own environment and that of 
the TLS to its own environment are weak enough that 
a Markovian master equation approach provides a good 
description of the dynamics. The combined qubit+TLS 
system has four quantum states. The quantity that we 
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Rabcd — 



9a(t) = I dLuS a (u)e-^\ (9) 

J — oo 

where S a (u>) is the background- noise power spectrum. In 
calculating R a bcd we neglect the imaginary parts, which 
renormalize the energy splittings of the qubit and TLS, 
and we assume that those corrections are already taken 
into account in our initial Hamiltonian. We do not use 
any secular approximation to simplify the tensor R a bcd 
any further. One of the main reasons for avoiding the sec- 
ular approximation is that we shall consider cases where 
the coupling strength between the qubit and the TLS is 
very small, which results in almost degenerate quantum 
states, a situation that cannot be treated using, for ex- 
ample, the form of the secular approximation given in 
Ref. |23. 

Once we solve Eq. {7J and find the dynamics of the 
combined system, we can trace out the TLS degree of 
freedom to find the dynamics of the reduced 2x2 den- 
sity matrix describing the qubit alone. From that dy- 
namics we can infer the effect of the TLS on the qubit 
decoherence and, whenever the decay can be fit well by 
exponential functions, extract the qubit dcphasing and 
relaxation rates. 



consider is therefore the 4x4 density matrix describing 
that combined system. We follow the standard procedure 
to write the Bloch-Redfield master equation as (see e.g. 
Ref. 0): 

Pab = —i^ab Pab + ^ Rabcd Pcd, (7) 
cd 

where the dummy indices a, b, c and d run over the four 
quantum states, io a b = (E a — Eb)/h, Ei is the energy 
of the quantum state labelled by i, and the coefficients 
Rabcd are given by: 



(8) 



IV. ANALYTIC RESULTS FOR THE 
WEAK-COUPLING LIMIT 



We first consider a case that can be treated ana- 
lytically, namely that of a strongly-dissipative weakly- 
coupled TLS. That is exactly the case where the tradi- 
tional weak-coupling approximation is expected to work. 
Here we perform a perturbative calculation on Eq. Q 
where the coupling strength A is treated as a small pa- 
rameter in comparison with the decoherence times in the 
problem. We shall discuss the differences between the 
predictions of the two approaches in this section, and we 
shall show in Sec. V that our results have a wider range 
of validity than the traditional weak-coupling approxi- 
mation. 

In the first calculation of this section, we consider the 
zero-temperature case. If we take the limit A — > and 
look for exponentially decaying solutions of Eq. J7J with 
rates that approach the unperturbed relaxation and de- 
phasing rates r[ q ' and we find the following approx- 
imate expressions for the leading-order corrections: 



a=q,TLS I In 



+9a(-t) 



5 a c^{d\&i a) \n){n\aj i a) \b)e 1Unit + (a\^ z a) \c)(d\ai a) \b)e 1 ' 
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which are expected to apply very well when r 2 TLS '' + 3> F^. The above expressions can be compared with those 
given in Ref. |T^ : 



(TLS) 



r 2 TLS) +(E q -E TLS ) 2 



di 2 - 2 di l+ p(TLS) SCCh ^2fc fl rJ' 



The two approaches agree in the limit where they are 
both expected to apply very well, namely when the de- 
coherence times of the TLS are much shorter than those 
of the qubit [note that we are taking T = in Eq. I|10|) ]. 
Our results can therefore be considered a generalization 
of those of the traditional weak-coupling approximation. 
We shall discuss the range of validity of our results in 
Sec. V.C.I. 

We now turn to the finite temperature case. In ad- 
dition to treating A as a small parameter, one can also 
perform a perturbative calculation to obtain the tem- 



perature dependence of the decoherence rates in the low 
temperature limit. In the general case where the qubit 
and TLS energy splittings are different and no assump- 
tion is made about the relation between qubit and TLS 
decoherence rates, the algebra is rather complicated, and 
the resulting expressions contain a large number of terms. 
Therefore we only present the results in the case where 
E q = Etls = E, which is the case that we shall focus on 
in Sec. V. In that case we find the additional corrections 
to the relaxation and dephasing rates to be given by: 



ip(q) _ 
oi 1T — 



oT (q) = \ 2 e~ E/kBT ( 4 cos2 flq cos2 ^tls sin 2 fl q sin 2 9 t -ls^i 



p(TLS) 
1 1 




(12) 



Note that the first term inside the parentheses agrees 
with the expression given in Eq. I jllj l for the TLS con- 
tribution to the dephasing rate. If the decoherence times 
of the TLS are much shorter than those of the qubit, 
the second term is negligible. A similar situation occurs 
when E q ^ Stls, i-e., all the terms can be neglected 
except for the one given in Eq. (fill) . 



V. NUMERICAL SOLUTION OF THE MASTER 
EQUATION 

We now turn to the task of numerically analyzing the 
effect of the TLS on the qubit with various choices of pa- 
rameters. We analyze any given case by first solving Eq. 
(|7|) to find the density matrix of the combined qubit+TLS 
system as a function of time. We then trace out the TLS 
degree of freedom to find the (time-dependent) density 
matrix of the qubit alone, which is perhaps most easily 
visualized as a curve in the Bloch sphere |2q. We then 
use the Hamiltonian of the qubit including the mean- 
field correction contributed by the TLS as a pole of ref- 
erence in the Bloch sphere, from which we can extract 
the dephasing and relaxation dynamics of the qubit. In 
other words, we transform the qubit density matrix into 



the qubit energy eigenbasis, such that the diagonal ma- 
trix elements describe relaxation dynamics and the off- 
diagonal matrix elements describe dephasing dynamics. 
The relaxation rate is then defined as the rate of change 
of the diagonal matrix elements divided by their distance 
from the equilibrium value. The dephasing rate is defined 
similarly using the off-diagonal matrix elements |2{j . 

Since our main goal is to analyze the different possi- 
ble types of behaviour in the qubit dynamics, we have 
to identify the relevant parameters that determine the 
different behaviour regimes. As discussed in Sec. II, 
the energy scales in the problem are the qubit and TLS 
energy splittings, their background decoherence rates 
(which are related to the environment noise power spec- 
trum), the qubit-TLS coupling strength and tempera- 
ture. Note that if the difference between the two en- 
ergy splittings is substantially larger than the coupling 
strength, the effect of the TLS on the qubit dynamics 
diminishes rapidly. The above statement is particularly 
true regarding the relaxation dynamics. We therefore 
consider only the case where the two energy splittings 
are equal, i.e. E q = -Etls- In other words, we take the 
TLS to be on resonance with the qubit. Furthermore, 
we take the energy splitting, which is the largest energy 
scale in the problem, to be much larger than all other 
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energy scales, such that its exact value does not affect 
any of our results. We take the temperature to be much 
smaller than E q , so that environment-assisted excitation 
processes can be neglected. We are therefore left with the 
background decoherence rates and the coupling strength 
as free parameters that we can vary in order to study the 
different possible types of behaviour in the qubit dynam- 
ics. 



A. Weak-coupling regimes 

Although the discussion of the criterion that distin- 
guishes between the weak and strong-coupling regimes is 
deferred to Sec. VI, we separate the results of this sec- 
tion according to that criterion. We start with the weak- 
coupling regimes. In Figs. 1 and 2, we show, respec- 
tively, the relaxation and dephasing rates of the qubit as 
functions of time for three different sets of parameters 
differing by the relation between the qubit and TLS de- 
coherence rates, maintaining the relation = 2r[ a \ 
where and are the background relaxation and 



dephasing rates, i.e. those obtained in the case A = 0, 
and the index a refers to qubit and TLS. The case A = 
is trivial, and we only show it as a point of reference to 
demonstrate the changes that occur in the case A ^ 0. 
All the curves shown in Figs. 1 and 2 agree very well 
with the formulae that will be given below. 



1. Relaxation dynamics 

Characterizing the dynamics is most easily done by 
considering the relaxation dynamics. From Fig. 1 we can 
see that in the case A 7^ 0, there are several possible types 
of behaviour of the qubit depending on the choice of the 
different parameters in the problem. As a general simple 
rule, which is inspired by Fig. 1(a), we find that for small 
values of A the relaxation rate starts at its background 
value and follows an exponential decay function with a 

characteristic time given by (r 2 TLS ' ) +r 2 q ^ — Ij after 
which it saturates at a steady-state value given by Eq. 
JTDJ, with E q = E TLS : 



dP cx (t)/dt 

P CX (t) ~ Pox (OO) 



-(q) 



A 2 sin 2 q sin 2 6- 



TLS 



2 r 



.(TLS) 



r(q) pv 
1 2 1 1 



(q)\ 



(i-ex P {-( r ™ + r 2 ^r^} 



(13) 



We can therefore say that the qubit relaxation starts with 
an exponential-times-Gaussian decay function for a cer- 
tain period of time, after which it is well described by an 
exponential decay function with a rate that incorporates 
the effects of the TLS, namely that given in Eq. (|10fl . 
Clearly the above picture is only valid when the expres- 
sion for the transient time given above is much smaller 
than l/I^ . Furthermore, Eq. (|T3|l is not well defined 



I2 — r^ qj = 0. However, even when the 



when r 2 TLS '' 

above condition about the short transient time is not sat- 
isfied, and even when the exponent in Eq. (|13[) becomes 
positive, we find that for times of the order of the qubit 
relaxation time, a very good approximation for the re- 
laxation rate is still given by Eq. i|13[l . The reason why 
that is the case can be seen from the expansion of Eq. 
(|13|l in powers of t: 

dput)/dt „ _ r ( q) _ 1 A2 sin2 0q sin2 W) (14) 



Pex(t) " Pcx(oo) 

which can be integrated to give: 



P ex (t) «exp{-rj q) i- A 2 sin 2 q sin 2 T Lsi 2 /4}, (15) 

where we have assumed that P e x(0) = 1 and that P ex (oo) 
is negligibly small. Eq. I|15f) describes the initial decay for 



any ratio of qubit and TLS decoherence times. Whether 
that function holds for all relevant times or it turns into 
an exponential-decay function depends on the relation 
between 1^ and r 2 TLS ^ + as discussed above. In 
particular, in the case when the TLS decoherence rates 
are much smaller than those of the qubit, Eq. H15|l holds 
at all relevant times, and the contribution of the TLS 
to the qubit relaxation dynamics is therefore a Gaussian 
decay function. We also note here that the relaxation 
rate shows small oscillations around the functions that 
we have given above. However, those oscillations have a 
negligible effect when the rate is integrated to find the 
function P ex (t). 



2. Dephasing dynamics 

The dephasing dynamics was somewhat more difficult 
to analyze. The dephasing rate generally showed oscilla- 
tions with frequency E q , and the amplitude of the oscil- 
lations grew with time, making it difficult to extract the 
dynamics directly from the raw data for the dephasing 
rate. However, when we plotted the averaged dephasing 
rate over one or two oscillation periods, as was done in 
generating Fig. 2, the curves became much smoother, 
and we were able to fit those curves with the following 
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FIG. 1: Relative corrections to qubit relaxation rate as a function of scaled time in the case of (a) strongly, (b) moderately 



and (c) weakly-dissipative TLS. The ratio r^ TLS) /rJ q) is 10 in (a), 1.5 in (b) and 0.1 in (c). The solid, dashed, dotted and 



dash-dotted lines correspond to X/V^ — 0, 0.3, 0.6 and 0.9, respectively. 8 q = tt/3 and #tls = Sn/8. 
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FIG. 2: Relative corrections to qubit averaged dephasing rate, i.e. after eliminating fast oscillations around the slowly-varying 
function, as a function of scaled time in the case of (a) strongly, (b) moderately and (c) weakly-dissipative TLS. The ratio 

The solid, dashed, dotted and dash-dotted lines correspond to X/V^ — 0, 
7r/3 and #tls = Stt/8. 



r (TLS)^ r (q) ig 1Q - n ^ Lg . n ^ and Q l . n ^ 
0.3, 0.6 and 0.9, respectively. 



simple analytic formula, which we obtained in an analo- 
gous manner to Eq. JT3J: 



1 

poi 



dp. 



01 



dt 



,(q) 



A 2 sin 6„ sin 



7TLS 



4 r. 



(TLS) 



i-ex P {-(ri TLS )-rW)t} 



(16) 



Starting from this point, the analysis of the dephasing 
dynamics is similar to that of relaxation. When the de- 
coherence rates of the TLS are much larger than those of 
the qubit, the dephasing rate starts from its background 
value but quickly reaches its steady-state value given by 
Eq. I|10|) . In the opposite limit, where the TLS decoher- 
ence rates are much smaller than those of the qubit, a 
good approximation is obtained by expanding Eq. 11161) 
to first order in t. In that case we find that: 



Poi(i) « Poi(0) exp [-T^t - A 2 sin 2 6 q sin 2 6W7s}- 

(17) 



B. Strong-coupling regime 

In the strong-coupling regime corresponding to large 
values of A, the qubit relaxation and dephasing rates 
as plotted similarly to Figs. 1 and 2 show oscillations 
throughout the period where the qubit is far enough from 
its thermal equilibrium state. Therefore one cannot sim- 
ply speak of a TLS contribution to qubit decoherence. 
Analytic expressions can be straightforwardly derived for 
the dynamics in the limit where the coupling strength is 
much larger than the decoherence rates. However, the 
algebra is quite cumbersome, and the results are rather 
uninspiring. Therefore, we shall not present such expres- 
sions here. 
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C. Further considerations 

1. Comparison with traditional weak-coupling 
approximation 

To demonstrate the differences between our results and 
those of the traditional weak-coupling approximation, we 
plot in Fig. 3 the relaxation rate (at the end of the tran- 
sient time) as a function of the coupling strength A. A 
similar figure can be obtained for the dephasing rate, but 
we do not include it here. Our numerical results agree 
with the perturbation calculation of Sec. IV (Eq. I1U|) 
up to the point where the coupling can be classified as 
strong, as will be explained in Sec. VI. Therefore we 
conclude that the results of our perturbation calculation 
have a much wider range of validity than those of the 
traditional weak-coupling approximation. The more non- 
negligible the TLS decoherence times are relative to those 
of the qubit, the larger the difference between the two 
approaches. Note that in our perturbation-theory calcu- 
lation we took the limit where A is much smaller than all 
the decoherence rates in the problem. It turns out, how- 
ever, that the results of that calculation are valid as long 
as A is substantially smaller than the TLS decoherence 
rates, assuming those are substantially larger than the 
qubit decoherence rates. No specific relation is required 
between A and the qubit decoherence rates in that case. 




FIG. 3: The steady-state qubit relaxation rate as a func- 
tion of the qubit-TLS coupling strength A. The solid line 
represents the numerical results, the dashed line is the 
perturbation-calculation result of Sec. IV and the dotted 
line is the result of the traditional weak-coupling approxi- 
mation given by Eq. I1H . 8 q — 37r/20, #tls = 27r/5, 

if VifLS) = .25, r< q) /rr LS) = i, and r< TLS Vr( TLS) = 2. 

For further demonstration of the differences between 
the predictions of the two approaches, we ran simulations 
of an experiment where one would sweep the qubit energy 
splitting and measure the relaxation and dephasing rates. 



We used a TLS with the parameters of Fig. 3(a) and A 
ranging from 0.15 to 0.5. In such an experiment, one 
would see a peak in the relaxation and dephasing rates 
at the TLS energy splitting. According to the traditional 
weak-coupling approximation, the height of the dephas- 
ing peak should be half of that of the relaxation peak. In 
the numerical simulations with the above parameters, we 
see a deviation from that prediction by about 25%. The 
relation between the shapes of the two peaks agrees very 
well with the expressions in Eq. (|10fl . That difference 
would, in principle, be measurable experimentally. Note, 
however, that since we are dealing with an uncontrollable 
environment, there is no guarantee that a TLS with the 
appropriate parameters will be found in the small num- 
ber of qubit samples available at a given laboratory. 



2. Two TLSs 

In order to establish that our results are not particular 
to single quantum TLSs, we also considered the case of 
two TLSs that are both weakly coupled to the qubit. If 
we take the two TLS energy splittings to be larger than 
the widths of their frequency-domain correlation func- 
tions, we find that the relaxation dynamics is affected by 
at most one TLS, depending on the qubit energy split- 
ting. The TLS contributions to pure dephasing dynam- 
ics, i.e. that unrelated to relaxation, are additive, since 
that rate depends on the zero-frequency noise. We then 
considered two TLSs with energy splittings equal to that 
of the qubit. We found that the TLS contributions to the 
qubit relaxation and dephasing dynamics are additive in 
both the large and small r( q ) /r^ TLS ^ limits. These results 
are in agreement with those found in Ref. |30|. where a 
related problem was treated. 



3. Entanglement 

It is worth taking a moment here to discuss the ques- 
tion of entanglement between the system and environ- 
ment. It is commonly said that in a Markovian master 
equation approach the entanglement between a system 
and its surrounding environment is neglected, a state- 
ment that can be misinterpreted rather easily. In order 
to address that point, we consider the following situa- 
tion: we take the parameters to be in the weak-coupling 
regime, where such a discussion is meaningful. We take 
the qubit to be initially in its excited state, with no en- 
tanglement between the qubit and the TLS. We find that 
the off-diagonal matrix elements of the combined system 
density matrix describing coherence between the states 
ItqJ-TLs) and UqtTLs) start from zero at t = and reach 
a steady state at the end of the transient time. Beyond 
that point in time, they decay with the same rate as the 
excited state population. We therefore conclude that the 
final relaxation rate that we obtain takes into account 
the effects of entanglement between qubit and TLS, even 
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though the density matrix of the qubit alone exhibits ex- 
ponential decay behaviour. 

VI. CRITERIA FOR STRONG COUPLING 
BETWEEN QUBIT AND TWO-LEVEL SYSTEM 

There are a number of possible ways one can define 
the criteria distinguishing between the weak and strong 
coupling regimes. For example, one can define a strongly- 
coupled TLS as being one that contributes a decoherence 
rate substantially different from that given by some weak- 
coupling analytic expression. One could also define a 
strongly-coupled TLS as being one that causes visible 
oscillations in the qubit dynamics, i.e. one that causes 
the relaxation and dcphasing rates to change sign as time 
goes by. We shall use the criterion of visible deviations in 
the qubit dynamics from exponential decay as a measure 
of how strongly coupled a TLS is. 

Even with the above-mentioned criterion of visible de- 
viations in the qubit dynamics from exponential decay, 
one still has to specify what is meant by visible devia- 
tions, e.g. maximum single-point deviation or average 
value of deviation. One also has to decide whether to 
use relaxation or dephasing dynamics in that definition. 
We have used a number of different combinations of the 
above and found qualitatively similar results. Those re- 
sults can essentially be summarized as follows: a given 
TLS can be considered to interact weakly with the qubit 
if the coupling strength A is smaller than the largest 
(background) decoherence rate in the problem. The ex- 
act location of the boundary, however, varies by up to 
an order of magnitude depending on which part of the 
dynamics we consider and how large a deviation from 
exponential decay we require. 

We have also checked the boundary beyond which our 
analytic expressions and numerical results disagree, and 
we found that the boundary is similar to the one given 
above. That result confirms the wide range of validity 
of our analytic expressions. Note in particular that even 
if the qubit-TLS coupling strength A is larger than the 



decoherence rates of the qubit, that TLS can still be con- 
sidered weakly coupled to the qubit, provided the TLS 
decoherence rates are larger than A. 



VII. CONCLUSION 

We have analyzed the problem of a qubit interacting 
with a quantum TLS in addition to its coupling to a 
background environment. We have characterized the ef- 
fect of the TLS on the qubit decoherence dynamics for 
weak and strong coupling, as well as weakly and strongly- 
dissipative TLSs. We have found analytic expressions 
for the contribution of a single TLS to the total deco- 
herence rates in the weak-coupling regimes, which is a 
much larger range than just the weak-coupling limits. 
We recover the results of the traditional weak-coupling 
approximation as a special case of our results, namely 
for a weakly-coupled strongly-dissipative TLS. We have 
found that weakly-coupled weakly-dissipative TLSs ex- 
hibit memory effects by contributing a non-exponential 
factor to the qubit decoherence dynamics. We have ver- 
ified that the contributions of two TLSs to the qubit re- 
laxation and dephasing rates are additive in the weak- 
coupling limit. We have discussed the transition from 
weak to strong coupling and numerically found that the 
transition occurs when the qubit-TLS coupling strength 
exceeds all the decoherence rates in the problem. 
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